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Abstract 

In this paper the nonperturbative analysis of the spectrum for one-particle excitations of the electron- 
positron field (EPF) is considered in the paper. A standard form of the quantum electrodynamics 
(QED) is used but the charge of the "bare" electron eo is supposed to be of a large value (eo ^> 1). 
It is shown that in this case the quasi-particle can be formed with a non-zero averaged value of the 
scalar component of the electromagnetic field (EMF). Self-consistent equations for the distribution 
of charge density in the "physical" electron (positron) are derived. A variational solution of these 
equations is obtained and it defines the finite renormalization of the charge and mass of the electron 
(positron). It is found that the coupling constant ao between EPF and EMF and mass mo of the 
"bare" electron can be connected with the observed values of the fine structure constant a and the 
mass of the "physical" electron m as follows (U = c = 1): 

1 m 
Qo ~ — ; mo~ — . 

a a 

It is also shown that although the non-renormalized QED corresponds to the strong coupling 
between EPF and EMF, the interaction between "physical" electron (positron) with EMF is defined 
by the observed value of the coupling constant a. It is proved that the translational motion of the 
"physical" particle is separated from its internal degrees of freedom. As a result the dependence 
of the one-particle excitation energy on its total momentum is defined by the formula E(P) = 
\/P 2 + m? and corresponds to the relativistic spectrum of a free particle with the observed value of 
mass m. Regularization of the terms of a series of the perturbation theory by a is due to the form- 
factor of the "physical" electron (positron) and corresponds to the cut-off momentum k max — mo. 
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I. INTRODUCTION 

Quantum electrodynamics (QED) is one of the basic elements of modern physics foundation. It allows to cal- 
culate the observed characteristics of fundamental electromagnetic processes with a unique accuracy. This means 
undoubtedly that the basic equations of QED correctly describe possible states of the electromagnetic (EMF) and 
electron-positron (EPF) fields. It is also well known that the problem of describing one-particle excitations of EPF 
can be formulated consistently in the framework of QED. The characteristics of these quasi-particles should determine 
the renormalization of the electron (positron) mass and charge because of their interaction with EMF. But the real 
solution of this problem in terms of the perturbation theory leads to some infinite values. An important property of 
QED is its renormalizability. It means that a special procedure can be developed which permits to exchange formally 
infinite parameters for their observed values when calculating the characteristics of electromagnetic processes. But 
at the same time Dirac wrote that "the calculation rules of QED are badly adjusted with the logical foundations of 
quantum mechanics and they can't be considered as the satisfactory solution of the difficulties" 1] (§81) . 

It's worth reminding that these rules are based either on some way of regularization for diverged integrals (for 
example, 0) or on adding counter-terms to Lagrangian of the system 8). As a result all observed characteristics 
are finite values and only unobserved values of the charge eo and mass m of the "bare" electron (positron) remain 
indeterminate because they include the diverged integrals. Recently a successive way for "finite QED" was formulated 
> [HI- I n this casc the perturbation theory series is calculated strictly for elements of the scattering matrix 
(S ~ matrix). The regularization of integrals is realized by means of accurately calculating the singular functions and 
the problem of describing the "physical" particle internal structure doesn't arise at all. 

Thus, from the "pragmatic" point of view the existing form of QED is a perfect technique for describing observed 
electromagnetic phenomena. Nevertheless, this approach is not completely accurate because "it is simply a way to 
sweep the difficulties under the rug" as Feynman wrote @]- 

What can stimulate searching for the solution of the problem of calculating the "bare" electron characteristics eo 
and mo in spite the fact that these values are unobservable and are not strictly included into experimentally measured 
physical values? First of all the possibility of their successive calculation is connected with the question whether the 
system of Maxwell and Dirac equations is closed and self-consistent when considering the processes involved with the 
interaction of electrons, positrons and photons? It is important to stress that we mean the internal consistency of QED 
only as the mathematical model for quantum field system. It can't be considered closed from the physical point of 
view because of its interaction with other quantum fields. In that sense QED in the existing form of renormalization is 
an unclosed theory because it includes an additional, external dimensional parameter which enables the regularization 
of integrals. 

Besides, the analysis of the "physical" electron structure could clarify the question about the continuity of space. 
The matter is that fundamental length was introduced in some papers as an additional parameter of the theory in 
order to realize the finite renormalization in QED. 

There are a lot of papers where the singularity-free electrodynamics is considered on the basis of a classical model 
for electron and vacuum polarization (see, for example, the more recent one |T^ | and references therein). But it is very 
important to consider the renormalization problem in the framework of the " logical principles of quantum mechanics" 
and relativistic invariance in order to develop some new nonperturbative approach in the field theory in application 
to the real physical system with completely defined Hamiltonian. At present nonperturbative methods are mainly 
studied for quite abstract quantum field models with a strong interaction (for example, [Tl)l. 

An essential role in motivating our work is played by the obvious contradiction between the small value of the 
coupling constant and its infinite calculated value showing a logical inconsistency of the perturbation theory as the 
main method used in QED. 

It is well known that the physical picture of the charge renormalization is mainly defined by the effect of polarization 
of electron-positron vacuum [l^. However, as the authors 01 (§128) noticed, it is "logically inconsistent" to analyze 
the connection between the charges eo and e in terms of relations based on the perturbation theory. So, in order to 
consider the renormalization problem we should reveal both the reason of divergence of the perturbation series and 
find the way for nonperturbative description of the system. 

In the present paper we suggest a possible variant for a nonperturbative investigation of the physical picture 
corresponding to the charge screening conditioned by the polarization of electron-positron vacuum. From the very 
beginning the initial charge is supposed to be larger than 1 ( eo ^> 1). It means that the nonrenormalized QED should 
be considered as the system of interacting fields with a str ong coupling. It was first shown by Pekar [15| and then 
analyzed more rigorously by Bogoluibov and Tyablikov [12], [23, that the main physical result of the strong coupling 
between the particle and the quantum field is defined by the self-localized (" polaronic" ) state of the particle. This 
state is due to the classical components selected from the quantum field because of its interaction with the particle. 
Later a lot of different forms of the quasi-particles arising as a result of a " polaronic" effect were described for many 
concrete physical systems. Another way for nonperturbative calculations in QED was considered recently in the paper 
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0. 

In this paper it is shown that the self-localized quasi-particle excitation of the EPF corresponds to the solution of 
QED equations for the case of eo ^> 1. The internal structure of such excitation is determined by a spatial charge 
distribution of EPF coupled with a classical scalar component of EMF separated by the same charge density. We 
associate such excitation with the "physical" electron (positron). It allows us to establish the analytical relation 
between QED initial parameters eo,mo with their observed values e,m and to justify the self-consistency of the 
original supposition. 

When describing the observed electromagnetic processes the relativistic invariance and renormalizability are the 
most important properties of QED confirmed by a lot of precision experiments. Therefore the self-consistent dynamic 
form- factors, which we have introduced for the internal structure of the "physical" electron, should be relativistically 
invariant and reveal themselves in the observed physical characteristics only in the form of the renormalized charge 
and mass of the electron at least up to very high energies of interacting particles. The analysis of this problem is 
considered in the second part of the paper. Our purpose is to substantiate that the supposition about the strong 
coupling between EPF and EMF doesn't logically contradict the standard perturbation theory in the renormalized 
QED. 

The paper is organized as follows. In Section 2 we consider the qualitative analysis of the physical picture of the 
electron-positron vacuum polarization when the quasi-particle excitation is formed. The use of quite simple scale 
estimation permitted us to find the characteristic correlations between the observed charge and mass of the electron 
and the initial parameters of the theory. Besides, the mathematical reason for the divergence of the perturbation 
series in QED is discussed. 

In Section 3 the method for nonperturbative calculation of the characteristics of self-localized one-particle excitation 
of the electron-positron field is developed in the framework of QED. The system of self-consistent equations for 
describing the charge distribution and the potential of the electrostatic field coupled with it is derived from the zero 
approximation of the approach under consideration. The total momentum of the excitation is supposed to be zero 
in this case. The variational solution of these equations is found and more precise relations between the charge and 
mass of " bare" electron and the corresponding values for " physical" electron are calculated. 

The dependence of the quasi-particle excitation energy E{P) on its total momentum is found in Section 4. As is 
known such dependence is rather complicated in the "polaron" problem and it is very difficult to find analytically 
its general form because the variables describing the internal dynamics of the system are not separated from the 
tranlational variables |2j|. Therefore the calculation of E{P) for QED is a very important verification for the non- 
perturbative approach in question when describing the observed properties of the "physical" electron (positron). It is 
shown in this section that nonperturbative QED satisfies all necessary conditions: the internal degrees of freedom are 
separated from the translational motion of the "physical" electron (positron) and its dispersion law E(P) corresponds 
to the standard dynamics of the relativistic free particle. 

One of the main purposes of our paper is differs greatly from most of the theoretical works where some results for 
experimental observation are predicted. On the contrary, we are going to prove that the internal structure of the 
"physical" electron considered in our paper doesn't reveal in its observable characteristics. Therefore, in Section 5 
the rules for transition from the initial theory with the strong coupling to the renormalized QED are justified. It is 
shown that the interaction between the " physical" electron and the transversal electromagnetic field is defined by the 
observed electron charge. It means that the characteristics of real electromagnetic processes can be calculated using 
the standard perturbation theory, regularization of the diverged integrals in QED being defined by the form-factor of 
the charge distribution in the "physical" electron (positron). 

II. QUALITATIVE ANALYSIS 

Let us first discuss the qualitative character of the processes which can determine the formation of the "physical" 
electron (PE). We use Dirac's interpretation for the ground state of the vacuum electron-positron field (VEPF) pj. 
As is known (for example, [LUp. such representation for the solutions of Dirac equation with negative energies is 
not rigorous enough but allows one to interpret intuitively the qualitative picture of the effects associated with the 
creation of pairs and polarization of VEPF in strong electromagnetic fields. 

It should be noted, that in this picture the ground state of VEPF is taken to be the energy zone completely filled 
by the "bare" electrons (BE) with the charge (— eo), eo > up to the negative maximal energy (— m ). The empty 
states of the continuous spectrum are separated from the border of the vacuum zone by the energy gap Ao = 2mo, 
with mo being the mass of BE. 

We will proceed from the above mentioned statement that QED is the closed model for two interacting fields in the 
sense that all the problems of renormalization in this theory are considered to be the result of some mathematical 
difficulties. It means that they should be overcome without including the interaction with other fields and without 
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any changes in the QED Hamiltonian for regularization. Thus, all further constructions are based on the assumption 
that the values m Q and eo are the only parameters of the theory, provided we use the system of units with K = c = 1. 

If the Dirac representation for the energy structure of VEPF is used, the one-particle excitation can be regarded as 
the appearance in the free energy zone of the "bare" electron localized near some point of space. Generally speaking, 
the assumption about local excitation breaks the translational symmetry of the system. It will be shown below 
(Section 4) that taking into account this symmetry one can find the dependence of the excitation energy on the total 
momentum and this energy spectrum proves to be relativistically invariant. Actually, in our qualitative consideration 
we discuss the quasi-classical picture for the formation of the state corresponding to " physical" electron (positron) at 
rest. To some extent this analysis is analogous to the estimation of the binding energy of the self- localized state of 
the electron in the ionic crystal which was considered in the first Pekar's work on the "polaron" theory [l5|. 

The electrostatic interaction between the excited BE and the electrons from a lower zone determines the polarization 
of VEPF which leads to some redistribution of the electrons in the filled zone and localized near the excited BE. In 
fact, it means that a surplus positive charge appears in this local domain and this charge can be described by some 
distribution of the "bare" positrons (BP). Neglecting the spin effects, the interaction between different charges in this 
stationary excited state can be defined mainly by the local self-consistent potential (LSP) of attraction ip(r) which 
corresponds to the classical part of the scalar component of EMF. It is supposed that the coordinate r is counted off 
from the center of the mass of the quasi-particle. 

Let the amplitude of LSP reaches its maximal value ip m in some point. Then the energy levels for the states 
of BE and BP near this point will be shifted down and up respectively along the energy scale by the values being 
proportional to the values eo|y m |- As QED is supposed to be the closed theory these shifts should be of the following 
form, taking into account following from the dimensional arguments: 



A±=mo/±(eo), (1) 
and the functions f±(&o) should be of the same form because of the charge symmetry. 

As a result, the real value for the energy gap between the energy of excitation and the top of the VEPF zone is 



A = A+ - A_ 

and it depends on the binding energy of the excited BE with the cloud of BP appearing due to the polarization of 
VEPF. 

The binding energy depending on the initial charge eo, it is possible to simulate the qualitative behavior of the 
dimensionless functions f± by the curves represented by the dotted lines in Figure 1. Let's denote by e$ that value 
of the BE charge at which the energy level of the one-particle excitation intersects locally the filled zone of VEPF. 
Evidently, the real value of eo is unknown a priori and we consider two "scenarios" of the renormalization which are 
qualitatively different. 

In the case of eo < ej, the sequence of the excited energy levels in the real system adiabatically corresponds to the 
structure of the energy zones in the ideal system. As a rule, it is one of the necessary conditions for the convergence 
of the series of the perturbation theory by the coupling constant Q . 

However, as the estimation given below shows, the real QED is apparently defined by another variant, i.e. eo > £q, 
corresponding to the strong coupling limit. In this case the quasi- intersection of the one-particle excited levels with the 
filled VEPF zone arises. As a result a new structure of the energy zones is to be formed upon including the interactions 
for "physical" particles. Fig.l shows this scheme with solid lines. An analogous scheme of the level "reconstruction" 
is considered for the effect of the pair creation in the homogeneous electric field 4] or in the Coulomb field of heavy 
atoms Similar behavior of the energy levels is observed in solid state physics when describing the "resonant" zones 
(for example, |6j). So, in the case of eo > ej one should expect the non-adiabatic correspondence between the levels 
of the " physical" electron-positron field and the basic set of the system levels without any interaction. 

The qualitative analysis of the possible behavior of the excited energy levels in QED allows us to establish two 
essential circumstances. Firstly, if the real value of eo is indeed in the range eo > e^, then the perturbation series 
for describing FEL structure, based on the set of states of the ideal system in the absence of interaction should be 
a diverged one because of the non-adiabatic change of the states when taking into account the interaction between 
EPF and EMF. 

Secondly, as Fig.l shows, the state corresponding to the top level of the filled (that is positronic) zone for "physical" 
EPF is actually the continuation of the bottom level of the empty (that is electronic) zone of VEPF if the level 
classification for non-interacting fields is used. It means that the quasi-particle excitation of VEPF in the "bare" 
electron form results in the state that should be considered as the "physical" positron (!). This statement has also 
a qualitative explanation. When the zones are locally overlapped due to the VEPF polarization, some part of BE 
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from the lower zone passes to free levels of the upper zone. Then they extend to infinity because of the repulsion 
with the initially negative fluctuation of VEPF. As the result, the screening cloud with the surplus positive charge 
is formed near the initial BE. Certainly, this intuitive arguments will be verified by analytical calculations in the 
following sections. 

Here we will consider a very simple and approximate quantitative analysis of the above qualitative picture of the 
renormalization. It is evident that the accurate description of the quasi-particle excitation of VEPF requires the 
solution of the related system of Dirac's equations for the electron-positron field and Maxwell's equations for the 
electromagnetic field (for example, 0) that will be considered below (§3). At the same time it is known that in the 
systems with Coulomb interaction there are no other dimensional values except the electron mass. In many cases 
it allows to find quite simple scaling estimations which define the physical characteristics of the system with the 
accuracy up to some numerical coefficients. So, in this section we use the quasi-classical scale analysis in order to 
define characteristic features of the "physical" electron or positron depending on the initial parameters of QED, taking 
into account the uncertainty relation but without considering the spin characteristics. 

In order to clarify the essence of this approach let's use it for estimating the energy of the ground state of the 
relativistic hydrogen-like atom with the nucleus charge Z. A similar estimation for a non-relativistic atom is well 
known, but it seems that it has not been considered for the relativistic case yet. If we consider the normalized wave 
function the quasi-classical formula for the relativistic energy of the hydrogen-like atom can be written as 



E H ~ V(M 2 +™ 2 , (2) 

a 

where the value a is the characteristic linear size of the domain where the electron is localized. It is related to the 
characteristic fluctuation of its momentum Sp by the uncertainty relation 



Here the fine structure constant is introduced 



u — Sp ~ — . (3) 
a 
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Taking the value u as the variational parameter it is possible to find the equation for it from the energy minimization. 
As the result the following estimation for the energy can be found 



u Zma 

= — Za = u — — =; 

V« 2 + m 2 y/1 - (Za) 2 

E H = mVl-H 2 . (4) 

It turns out that the quasi-classical estimation leads to the exact value for the ground state energy of the relativistic 
hydrogen atom as it does in the non-relativistic case. This coincidence may be accidental, but it should be remembered 
that it usually testifies the existence of closed classical trajectories in the system. The most essential point for us 
is that this scale estimation correctly defines that critical value of the nucleus charge at which Dirac energy gap 
disappears and the creation of the electron-positron pairs becomes possible [5| . 

Let us apply a similar estimation to considering the problem of the formation of the quasi-particle excitation. In 
this case an additional degree of freedom characterizes the physical structure of the quasi-particle together with the 
values of u and a. So we should introduce one more parameter for it. In fact, if the initial fluctuation of VEPF 
corresponds to BE its integral charge should be equal to the value (— eo). At the same time the integral screening 
charge q+ induced by the vacuum polarization may not be coincide with eo for a general case. The difference between 
these values defines the observed quasi-particle charge e which corresponds to "physical" electron. Thus wc introduce 
the additional variational dimensionless parameter x as the value defining the probability P + of finding the "bare" 
positron in the polarization cloud around the "bare" electron. As we consider the excitation normalized by one 
particle, then by definition 



p+ = t— ; P- = — T , (5) 

1 + 35 X + 1 



that is the value x = corresponds to the state with pure BE. 
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Then the charge balance in the quasi-particle excitation is 

x x — 1 

Q+ = t-; — e ; e=-(q--q + ) = — — -en. (6) 
1 + x x + 1 

Note that we don't fix the sign of the observed charge e relative to the charge of the excited BE (- eo). 

It is evident, that in the case of the standard perturbation theory the one-electron excitation of EPF is defined by 
the parameter x — 0. Let us indicate, however, that even the simple scale estimation shows that the minimal energy 
of the one-particle excitation is achieved with a non-trivial value of the parameter x. 

Suppose that the momentum fluctuations are equal for the electron and positron states forming the self-localized 
excitation. Then the quasi-classical estimation for the energy of the "physical" excitation can be written as follows if 
one takes into account the physical meaning of the parameter x 



u 2 + m 2 J—^- -ej—^- < <p(r) > -E = 0. (7) 
1 + x 1 + x 

Now we should estimate the average value of the self-consistent potential by means of the introduced parameters. 
The potential itself depends on the local distribution of the positive and negative charge densities ±enp±(r) in the 
quasi-particle and is defined by the Poisson equation which has the following form in the chosen system of units 

If no additional parameters for the local charge distribution are used, then the scale consideration leads to the 
following integral estimation for the averaged potential 

, % en fx — 1) en (x — 1) 
< ^ > ^t«i^a = -«(xTT ) U - ® 

Substituting the estimation J^J into Equation J2J we can find the closed variational formula for determining the 
unknown parameters x, u and energy E 



^Ju 2 + m 2 {l - x 2 ) + a (l - xfu - E{1 + xf = 0; 



a = ±. (10) 

Varying on the parameters leads to 

; " 2 (1 " A + "o(l - ^) 2 = 0. (11) 

- 2xJv? +m 2 + 2a Q (x - l)u - 2E(l + x) = 0. (12) 

The solution of Equation pi(l exists only if x > 1 and corresponds to the above given renormalization variant when 
eo > eg, so one can find 

moaot x — 1 



U v/1^2 

Substituting this solution into Eqs. lfT2*|l and ifTHIl leads to 



(13) 



2 m 
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a 



x 



2 t 2 =e\ I- alt 2 : 



-t + a^t 3 = eJ 1 - a^ 2 



o' 

E = m e, (14) 



with e being the dimensionless parameter defining the self-energy of the quasi-particle. 

Excluding e from the last two equations leads to the equation for the screening parameter x: 



2al(x- l) 2 = (x+ l) 2 . (15) 

This equation relates x to the initial charge eo which, in its turn, can be expressed by the observed charge e of the 
quasi-particle by means of the normalized condition 0. The explicit form for the values x and e is: 

t= — x = a ° t + 1 ; \fl- alt 2 = 
ao 



y/2' aoV2-l' V ° v/2 



E = -m t x /l-a 2 t 2 = -^. (16) 

zao 



Substituting these equations into the normalized condition © shows that according to the qualitative analysis the 
observed charge of the quasi-particle proves to be positive (e > 0), i.e. it corresponds to "physical" positron . This 
result is also consistent with the negative value of the self energy of the quasi-particle that defines the top of the 
" physical" positron zone (-m) . Its absolute value should be equal to the observed mass of the " physical" positron (see 
Fig.l). Thus, in the framework of the accuracy of the considered quasi-classical estimation the connection between 
renormalized and initial parameters of QED is defined by the formulas clearly indicating the nonperturbative character 
of the above calculation 



1 1 

a V2 2a 
m m 

m = - — ; m = — . (17) 
Zao a 

To describe the characteristics of the "physical" electron one has to repeat all the calculations considering "bare" 
positron as the initial fluctuation of VEPF, not changing the values of the renormalization parameters because of the 
charge symmetry. The numerical values of the initial parameters of QED are of some interest, too: 

a = ~ 68.5; 
la 

Tfl 

m = — ~ 70MeV. (18) 
a 



III. VARIATIONAL EQUATIONS FOR ONE-PARTICLE EXCITATION IN QED 

As mentioned above, the Dirac interpretation of the energy levels of the vacuum electron-positron field is quite 
qualitative and the estimation obtained on its basis in Section 2 is useful mainly for visualizing the formation of the 
quasi-particle excitation, corresponding to the "physical" particle. Therefore it is important to answer the question: 
whether there exists the configuration of the interacting fields which leads to the analogous structure of one-particle 
excitation in the framework of the accurate formulation of QED? 

Let's use the hamiltonian QED in the Coulomb gauge [14j, This representation is most convenient for us because 
it extracts the electrostatic field which, by definition, makes the main contribution to the vacuum polarization in the 
limit of strong coupling: 
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H = J dr{$*(r)[a(f+ e A(f}) + (3m }if>(r) + e £{r)p{r) - ~{V(p{r)) 2 } + ^u;(£)% A ; 

kX 

P(r) = \[^{f)i>(f) - $(r)r(r)}. (19) 

We suppose here that the field operators are given in the Schrodinger representation, the spinor components of the 
EPF operators being defined in the standard way 0] 

= E / J^^H u ^ e ~ tpP + b fs v*_ P _ su e^}. (20) 

In these formulas a, (3 are Dirac matrixes; up sv and vp SI/ are the components of the bispinors corresponding to the 
solutions of Dirac equations for the free "bare" electron and positron with the momentum p and spin s; ap s (at s ) and 
bps{bp s ) are the annihilation (creation) operators for the "bare" electrons and positrons in the corresponding states. 

The field operator A(r) and the operator of the photon number are related to the transversal electromagnetic 
field and their explicit form will be written below. It should also be reminded that the parameter eo was introduced 
as the positively defined value. 

We should especially discuss the dependence of the QED hamiltonian on the operator of the scalar field 



!fi(r) = j dky % e i%? . (21) 

It is known that these operators could be excluded from the hamiltonian in the Coulomb gauge. For that 
purpose one should use the solution of the operator equations of motion for (p? supposing that the "bare" electrons 
are point- like particles and referring to "self-action" is equivalent to the substitution of the initial mass for the 
renormalized one. In the result the terms with scalar fields in the hamiltonian are reduced to the Coulomb interaction 
between the charged particles. But we can't use this transformation of the hamiltonian i|19|) because it is the dynamics 
of mass renormalization that we investigate. 

One more problem is connected with a negative sign of the term corresponding to self-energy of the scalar field. 
If the non-relativistic problems were considered then the operator of the particle kinetic energy would be positively 
defined and the negative operator with the square-law dependence on (p(r) would lead to the "fall on the center" [3(| 
as the energy minimum would be reached at an infinitely large field amplitude. However, if the relativistic fermionic 
field is considered then the operator of the free particle energy (the first term in formula H19|) ~) is not positively 
definite. Besides, the states of the system with the negative energy are filled. Therefore, the stable state of the system 
corresponds to the energy extremum (!) (the minimum one for electron and the maximum one for positron excited 
states) . It can be reached at the finite value of field amplitude (see below) . The same reasons permit one to successfully 
use the states with indefinite metric 0) in QED although it leads to some difficulties in the non-relativistic quantum 
mechanics. 

According to our main supposition about the large value of the initial coupling constant eo we are to realize the 
nonperturbative description of the excited state which is nearest to the vacuum state of the system. The basic method 
for an nonperturbative estimation of the energy is the variational approach with some trial state vector |<i>i > for an 
approximate description of the one-particle excitation. The qualitative properties of the self-consistent excitation in 
the strong coupling limit [l5l | show that such trial vector should correspond to the general form of the wave packet 
formed by the one-particle excitations of the "bare" EPF. Besides, the effect of polarization and the appearance of 
the electrostatic field ip(r) should be taken into account so we consider |<!>i > to be the eigenvector for the operator of 
the scalar field. Now, let's introduce the following trial vector depending on the set of variational classical functions 
; <p{r) for an approximate description of the quasi-particle excited state of the system: 

|$! >~ \^\U q ^ s ^(r)) >= J dq{U q ~ s a+ s + V^ s b+ s }\0;0 ;l p(r) >; 
0(f)|O;O; >=<p(?)\0;0;<p(f) >; a 3 - s |0; 0; tp{f) >= 6 g - s |0; 0; ip(f) >= 0. (22) 
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The ground state of the system is |<j>o >= |0; 0; >, if we use the same notation. It corresponds to the vacuum of 
both EPF and EMF. 

Firstly, let's consider the excitation with the zero total momentum. Then the constructed trial vector should satisfy 
the normalized conditions resulting from the definition of the total momentum P and the observed charge e of the 
"physical" particle: 

< >= ]T dqq[\U<? a \ 2 + |^| 2 ] = P = 0; 

S 

£ dq[\U qs \ 2 + \V qs \ 2 ] = l. (23) 



< ^\Q\^ >= e ]T dg\\V qs \ 2 - \U qs \ 2 } = e. (24) 

s 

The condition 123|) restricts the functions U qs ; V qs in the sense that they depend on the modulus of the vector q 
only. It should also be taken into account that the trial vector > is not the accurate eigenvector of the exact 

integrals of motion Q and P as it represents the accurate eigenvector of the hamiltonian |$i > only approximately. 
Therefore, in the considered zero approximation the conservation laws for momentum and charge can be satisfied only 
on the average, and this leads to the above written normalized conditioned. Generally, the equation (|24|l should not 
be considered as the additional condition for the variational parameters but as the definition of the observed charge of 
the "physical" particle if the initial charge of the "bare" particle is considered as the given value. Therefore the sign 
of the observed charge is not fixed a priori,the same was the case for the qualitative estimation in iffijl. Calculating 
the sequential approximations to the exact state vector |$i > (see §4) should restore the accurate integral of motion 
as well. An analogous problem appeares in the "polaron" theory when the momentum conservation law was taken 
into account for the case of the strong coupling (for example, [22], pfij)- 

In this respect the variational approach differs essentially from the perturbation theory where the zero approximation 
for one-particle state is described by one of the following state vectors: 



|$i >~ |$i PT) e>=atj0;0;0>; |$ x >~ \<S>[ PT) p >= fetJO; 0; > . (25) 

These vectors don't depend on any parameters and are eigenvectors of the momentum and charge operators. But 
they correspond to one-particle excitations determined by the charge eo of the "bare" electron and the field ip(r) = 0. 
As supposed the introduction of the variational parameters into the wave function of the zero approximation enabled 
us to take into account the vacuum polarization, but the variational approach in this form allows to conserve the 
exact integral of motion only on the average. This problem will be considered in more detail below (Section 4). 

So, the following variational estimation for the energy Ei(P = 0) of the state corresponding to the "physical" 
quasi-particle excitation of the whole system is considered in the strong coupling zero approximation: 



E 1 (0) ~ E^[U qs ;V qs ;ip(r)} =< ^\H\^ >, (26) 

where the average is calculated with the full hamiltonian (|19f) and the functions U qs ; V qs ; ip(r) are to be found as 
the solutions of variational equations 



dE[°\u qa ;V qa -M^) = dE^_ = dE^ = 

dU qs dV qs d<p(f) [ ' 

under some additional conditions (|23I - I24|I . 

It is quite natural, that the ground state energy is calculated in the framework of the considered approximation as 
follows 

E ~ 4 0) =< *o|^|$o > • (28) 
Then the energy of the quasi-particle representing the "physical" electron (positron) is calculated by the formula 
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E(P = 0) = E{ 0) -4 0) - ( 29 ) 

Further discussion is needed in connection with the application of the variational principle (|26I27I) for estimating the 
energy of the excited st ate, because usually the variational principle is used for estimating the ground state energy. As 
far as we knew Caswell |l6j was the first to apply such method for nonperturbative calculation of the excited states of 
the anharmonic oscillator with an arbitrary value of the anharmonicity . He called this approach as " principle of the 
minimal sensitivity" . It was shown in our paper jl7j | , that the application of the variational principle to the excited 
states is actually the consequence of the fact that the set of eigenvalues for the full hamiltonian does't depend on 
the choice of the representation for eigenfunctions. As the result, the operator method (OM) for solving Schrodinger 
equation was developed as the regular procedure for calculating further corrections to the zero-order approximation. 
Later OM was applied to a number of real physical systems and proved to be very effective when calculating the 
energy spectrum in the wide range of the hamiltonian parameters and quantum numbers ( |l8j|. |l9j |. |2f)j | and the cited 
references) . These results can serve as the background for the application of the variational estimation to the excited 
state energy in QED. The procedure for calculating the successive approximations is considered in Section 5. 

The average value in Eq. i|26[l is calculated neglecting the classical components of the vector field (quantum 
fluctuations should be taken into account in the further approximations) 



< *f ) |^*(f)[3i!(» i )]^(T ! )|*i 0) >= 0. (30) 

It should be noted that the possibility of constructing self-consistently the renormalized QED at the non-zero 
vacuum value of the scalar field operator was considered before |27j but the solution itself of the corresponding 
equations was not discussed. 

It is also essential that the normal ordering of the operators || wasn't introduced in the hamiltonian l|19fl . Therefore, 
its average value includes the infinite values corresponding to the vacuum energies of EPF for both the ground and the 
excited states. However, if the normalized conditions (|23|l are fulfilled these values are compensated in the expression 
(|29|l for the quasi-particle energy. Then the functional for determining the zero approximation for the energy of the 
EPF one-particle excitation is defined as follows: 



E(P = 0) = dr 



dq 



(2tt) 3 / 



'2 



^ ' s,s' n,v 

e ip(r)S^]U qs u ?SL , 



~ f df[V^(f)] 2 . (31) 



In order to vary the introduced functional let us define the spinor wave functions (not operators) which describe 
the coordinate representation for the electron and positron wave packets in the state vector \<§>\ 



(0) >: 



e iqr . (32) 



In particular, if the trial state vector is chosen in one of the forms l|25[) of the standard perturbation theory, the 
wave functions (|32|) coincide with the plane wave solutions of the free Dirac equation. For a general case the variation 
of the functional i|31[l by the scalar field leads to 



E(0) : / dr{^ + (r)[{-iaV + f3m ) + ^e (p(r)}^ (r) - 
V +c (r)[(-iaV + pm ) + -e ^(f)]* c (f); 

he/,-?' WrC/-/ 



df[^+(r)^(r) + * +c (r')* c (r')] = 1; (33) 
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J df[^ + {f)^(f) + * +c (r')* c (r')] = 1. (34) 

The principal condition for the existence of the considered nonperturbative excitation in QED is defined by the 
minimum of the functional (|33|l corresponding to a non-zero classical field. The structure of this functional shows 
that such solutions of the variational equations could appear only if the trial state vector included the superposition 
of the electron and positron wave packets simultaneously. 

The functional (|33|l is the analogue of the quasi-classical expression © . Equation l|33(l and the Fourier representation 
1)22( 1 for the trial vector clearly indicate that the supposition concerning the localization of the functions ^f(r) near 
some point doesn't contradict the translational symmetry of the system because this point itself can be situated at any 
point of the full space with equal probability. The general analysis of the correlation between the local violation of the 
symmetry and the conservation of accurate integral of motion for the arbitrary quantum system was first discussed 
in detail by Bogoluibov in his widely known paper "On quasi-averages" |2lj . Analogous analysis in the considered 
problem we'll discuss in Section 4. 

Varying the functional ((33)1 by the wave functions *f?(r) $> c (f) taking into account their normalization leads to the 
Dirac equations in their usual form describing the electron (positron) motion in the field of potential f(r): 



{(-iaV + l3m Q ) + e ip(r) - E{0)}^(r) = 0; 

{(-iaV + (3m Q ) + e <p(r) + £(0)}* c (r) = 0. (35) 

We should discuss the procedure of separatiing variables in more detail, because of the non-linearity of the obtained 
system of equations for the wave functions and the self-consistent potential. Since the considered physical system has 
no preferred vectors if P — 0, it is natural to regard the self-consistent potential as spherically symmetrical. Then 
the variable separation for the Dirac equation is realized on the basis of the well known spherical spinors : 



Here fljiM are the spherical spinors describing the spin and angular dependence of the one-particle excitation 
wave functions; j, M are the total excitation momentum and its projection; the orbital momentum eigenvalues are 
connected by the correlation I + I' = 2j. It is natural to consider the state with the minimal energy as the most 
symmetrical one, corresponding to the values j — 1/2; M = ±1/2:1 = 0; I' = 1. This choice corresponds to the 
condition that in the non-relativistic limit the "large" component of the spinor $ ~ g corresponds to the EPF 
electronic zone . In this case the unknown radial functions /, g satisfy the following system of the equations: 



^y— - ~(rg) -(E + m - e ip(r))(rf) = 0; 
dr r 

^1 + -(rf) + (E- m - e^{r))(rg) = 0. (37) 
dr r 

The states with various projections of the total momentum should be equally populated in order to be consistent 
with the supposition of the potential spherical symmetry with the equation 1(34( 1 . So, the total wave function of the 
"electronic" component of the EPF quasi-particle excitation is chosen in the following form: 



xt = ^[V,i/2 + fi i/2,!,-i/2]; i = 0;l. (38) 



In its turn, the wave function is defined on the basis of the following bispinor: 



*W = 



-ifi(r)tljiM 



(39) 
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The radial wave functions f± , gi in this case satisfy the following system of equations 

+ -(rgi) -(E + m + <w(r))(r/ x ) = 0; 
dr r 

- -(r/x) + (E - m + e cp(r))(r 9l ) = 0. (40) 
dr r 

These equations correspond to the positronic branch of the EPF spectrum with the " large" component ~ g\ in the 
non-relativistic limit. 

It is important to note that the functions \P and W c satisfy the equations l|33|l with the same energy E(0). It imposes 
an additional condition of the orthogonality to them: 

<* c |*>=0. (41) 

Taking into account this condition and also the requirement that the states with different values of M should be 
equally populated one finds the "positronic" wave function 



1 \x5j c -Mi c ]-( -i /i( r )Xo 

^^1/2,0,1/2 W l/2,0,-l/ 2 J - ^ gi ^ x - 



Xj = ^[ni/2,i,i/ 2 - ^1/2,1,-1/2]; i = 0;l. (42) 

The equation for the self-consistent potential follows from the definition of <p(r) in formula (|35|) taking into account 
the normalization of the spherical spinors [7J: 

The boundary condition for the potential is equivalent to the normalization condition 124(1 and defines the charge 
of the " physical" electron (positron) e 

VWIr^oo = J— = 
47IT 

poo 

^ J o ^r 1 [/ 2 (r 1 ) + <? 2 (^i)-/i 2 (^i)-<??(n)]- (44) 

It is important to stress that the form of the functions given above is defined practically uniquely by the imposed 
conditions. At the same time the obtained equations are consistent with the symmetries defined by the physical 
properties of the system. The first one is quite evident and relates to the fact that the excitation energy doesn't 
depend on the choice of the quantization axis of the total angular momentum. 

Moreover, these equations satisfy the condition of the charge symmetry |7j. Indeed, one can check directly that one 
more pair of bispinors leads to the equations completely coinciding with (|37|l . (|40|) 



^ = V-/iV)4mJ ; (45) 



%m - ( ) . (46) 



jim y -ig{r)n jm 

It means that these bispinors allows one to find another pair of the wave functions which are orthogonal to each 
other and to the functions J3HJl,l|I2Jl but including the same set of the radial functions 

q, = ( igi(r)xl V ^ c = ( -f(r)xt \ m 



14 



These functions differ from the set (|38(l . 1421) in that they have another sign of the observed charge of the quasi- 
particle due to the boundary condition (|44p. 

Let us now proceed to the solution of the variational equations. It follows from the qualitative analysis that 
the important property of the trial state vector freedom is the relative contribution of the electronic and positronic 
components of the wave function. Therefore let us introduce the variational parameter C in the following way: 



r 2 dr[f 2 (r)+g 2 (r)} = -±-; 
o 1 + L 

C 



r z dr[ft{r)+ai{r)] = - 



10 x T C 

The dimensionless variables and new functions can be introduced 



(48) 



rm ; E = em ; e ip(r) = m <p(x); 



47T 



u(x)y/mu~ = rg{r); v(x)^/m^ = r/(r); Ux{x)^/m^ = rg x (r); v 1 (x) y /rn^ = r/i(r). 



(49) 



As the result the system of equations for describing the radial wave functions of the EPF one-particle excitation and 
the self-consistent potential of the vacuum polarization can be obtained. These equations can be clearly interpreted 
and the equations (|10I - I12(I serve as their quasi-classical analogue. 



du 

dx 
dv 

dx 



u — [e - 

x 

1 

+ -v+ (e- 



(f>{x))v 
(j}(x))u 



du-\ 1 , , . „ „ 

+ _ Ul - ( e + 1 + (j> (x)) Vl 
dx x 



dvi 
dx 

(j>(x) = a Q 



1 



vx + (e — 1 + cp(x))ui 



dy 



p(y) 



l 



dyp(y)} 



p(x) = [u 2 (x) + v 2 {x) - u\{x) - v\{x)\. 



(50) 



These equations should be complemented by the conditions which correspond to the relations between the observed 
charge and mass, by means of their initial values and eigenvalues e,C of the equations (|50|) : 



£(0) = emo = ±to; 
(1-C) 



eo 



(1 + C) 



= ±e. 



(51) 



The sign under the renormalization conditions (|51|) will depend on the correlation between the value eo, calculated 
on the basis of the equations l|5U|l and the critical parameter eg as it was discussed in Section 2. 

The mathematical structure of the non-linear equations (|50|l is in m any respects similar to that of the equations 
determining the ground state of the system in the "polaron" problem 15], but the number of the equations in our 
case is larger. Generally speaking, we can find the direct numerical solution of these equations (|23|). But in the 
present paper we are restricted by an approximate variational solution which is usually close to the exact solution in 
the similar problem |15| . In order to find this solution let us use the functional which is an exact consequence of the 
functional ^ for E(0): 



/>OC 

I = / dxUu'v - v'u) - 2 (e+ 1 - (j))v 2 - (e - 1 - (j))u 2 + 

Jo x 



(u[vi - v[ui) + 2^- - (e + 1 + cj))vl - (e - 1 + c/>)u 2 } = 0. 



x 



(52) 
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The potential </>(x) in the functional i|50[) is supposed to be calculated as the exact solution of the Poisson equation, 
as it was made by Pekar |2g • The choice of the form for trial functions is defined by the conditions of the analytical 
behavior for the solutions in the limit x — > and their asymptotic for x — > 00. The simplest set of the functions 
satisfying these conditions and normalized in accordance with (|48|l is as follows: 



u{x) = 2AJ 7(1 ^ c) (1 + bx)xe- bx ; 



v(x)=2bx\ ^- >" e~ bx ' 



m{x) = 2bx 2 J±— y~e 



r l ( y) =2A 1 J (1 ^ c) (l + M, f ''. ( •*>-"> ! 



A 1 ,- 


A 2 )b 3 


' 3(1 


+ C) ' 


\c- 


A\)b 3 


3(1 


+ C) 



-bx . 



Unlike the quasi-classical estimation (J7J) this set of functions depends on 4 variational parameters (A, A\ , C, b) with 
quite evident physical interpretation. The self-consistent potential is calculated by means of the last equation in l|5U|l : 



{(1 - C) e- 2ba; [42(l - C) + 3bx(2l - 2A Z + 2A\ - 21C) 



.t(1 + ) l ' 7 42 

66 2 x 2 (7 - 2A 2 + 2A\ - 7C) + 2b 3 x 3 (7 - 4A 2 + AA\ - 7C)]}. (54) 



As the result the functional Ij52(l transforms as follows 



I=(C-l)-e(C + 



1) + 2b^(A iy Jc-A 2 - A^l-A 2 ) + 2(A 2 - A\) 



In a general case the variation of the dimensionless energy e in equation 1(55(1 by 4 parameters leads to rather 
complicated equations. Their numerical solution depends on the value «o which, in its turn, should be calculated 
from the renormalization condition 151(1 . However, the observed charge of the "physical" electron is defined by the 
parameter qq > 1 (see the final result and the quasi-classical estimation). In this limit the unknown solution are 
simplified and can be represented in the analytical form 

A 2 -A\^, (l-O*-; 

e= gft*M , (56) 

where the introduced parameters are defined in terms of the conditions of numerical functional minimum q(t, z, A, b) 
which doesn't depend on the coupling constant ag for this approximation: 



/ 3 til — 2A 2 ) 
q(t, z,A,b) = 2t-z- bJ 7{1 _ A2) [-^ + Az] + 

,,9t 2 3tz 93z 2 N 

6( 784 + 56 + W } (5?) 

When 3 parameters are excluded from the corresponding variational equations the algebraic equation for the last 
one can be obtained 

1519 - 12376A 2 + 1156CL4 4 + 1088A 8 = 0. (58) 
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A numerical solution of this equation allows to calculate the other parameters and the value of the functional q = qa 
at the point of minimum 

A = -0.920843; 6= 1.48024; z = -1.12751; 

t = -0.956075; q Q = -0.784642. (59) 

It can be noted that the choice of a negative sign for A value follows from the condition that the parameter b should 
be positive. 

The signs of the parameters z, qo show that properties of the quasi-particle resulting from the one-electron excitation 
of EPF correspond to the "physical" positron as it was obtained for the quasi-classical analysis. But the numerical 
values of the initial QED parameters expressed through the observed values by means of Eqs. i|39|) slightly differ from 



ao = — ^ 174; a = — ~ — ; 
a 47r 137 

m = ^-m ~ 222to ~ 113(MeV); 
\lo\ 

r° e = — h = l -H-^0A16r e . (60) 

TOoO 0Z Z TO 

The last equation shows that the " electromagnetic radius" r® of the quasi-particle is of the same order of magnitude 
as the parameter for the "classical electron radius" r e in the renormalized QED. Besides, it is worth mentioning that 
the mass of "bare" electron seems to be close to the observed mass of muon (106 MeV). 



IV. EPF QUASI-PARTICLE EXCITATION WITH ARBITRARY MOMENTUM 

In the previous section we have considered the possibility for the resting quasi-particle with a non-trivial self- 
consistent charge distribution, the finite energy E(0) and a zero total momentum P = to exist in the framework of 
nonperturbative QED. 

The obtained solution allows one to imagine the internal structure of the resting "physical" electron (positron) as 
the strongly coupled state of charge distributions of the opposite sign, the large values of integral charges of these 
distributions compensating each other almost completely and their heavy masses are being " absorbed" by the binding 
energy. 

Actually the energy ±E(0) defines the boundaries of the renormalized electron and positron zones resulting from 
the strong polarization of EPF when the excitation appears (Figure 1). But this excitation could be really interpreted 
as the "physical" electron (positron) if the sequence of the levels in every zone determined by the vector (Fig.l) 
is described by the relativistic energy spectrum of real particles, that is 



E(P) = y/P 2 + E 2 (0) = ^P 2 +m 2 . (61) 

Only in this case the energy E(0) can really be used for calculating the mass too of the "bare" electron. 

It is worth saying, that the problem of investigating the dynamics of the self-localized excitation should be solved for 
any system with the strong interaction between quantum fields in order to calculate its effective mass. For example, 
a similar problem for Pekar "polaron" in the ionic crystal was considered in the papers j^, j^, [3], |2^|, pfij 
and in many more recent works. It is essential that because of the non-linear coupling between the particle and a 
self-consistent field the energy dispersion E{P) for the quasi-particle proved to be very complicated. As the result, its 
dynamics in the crystal is similar to the motion of the point "physical" particle only at rather small full momentum. 

However, in the case of QED the problem is formulated fundamentally in a different way. At present there is 
no doubt that the dynamics of the "physical" electron is described by the formula l|61|l for any (!) values of the 
momentum P. It means that the considered nonperturbative approach for describing the internal structure of the 
" physical" electron should lead to the energy dispersion law (16111 for the entire range of the momentum P. The analysis 
of this problem can be considered as an important test for finding out whether our method suits the experimental 
data. 

In order to describe the quasi-particle excitation with the non-zero momentum in the framework of the variational 
principle, let us introduce a more general trial state vector than the one in the formula (|22|l 
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>^ \^\u ?Ps ;V, Ps ;<p(r,P):P) > = 
/ dq{U w _ Pls a++V w _ Pls b+}\0;0^(r,P);P>. (62) 

We should stress once more that the wave packet ()62|) doesn't distinguish any point in the space when the energy 
of the stationary state is estimated by the variational functional. The matter is that the shift of the coordinate in the 
self-consistent potential ip(f, P) is reduced to the change of the variable when integrating in the full QED hamiltonian. 
The constructed state vector differs from the one-particle wave packet (|22H in that the momenta of the EPF excitations 
forming the internal structure of the " physical" electron are counted off from the vector P. The variational parameters 
£/_. p s ■ V~ p s depend only on the modulus of the vectors difference because of the normalization condition and the 
conservation on the average of the integral of motion 

< ^\P\^ >= £ dqq[\U w _ Pl f + \V w _ P{ f] = P; 

S 

£ dg\\U qs \ 2 + \V qs \ 2 ] = l. (63) 



< $( 0) |Q|$< 0) >= e ^ dq[\V qs \ 2 - \U qs \ 2 } = e. (64) 

s 

Then the analogue of the functional H31(l determining the energy E(P) — — Eq ^ of the quasi-particle excitation 
with an arbitrary momentum is defined as follows 

E{P) = J dfj J Y^ U \f-?W u 9>>tM& + P m °^» + 

e (p(r, P)S tlv ]U^_p ls u^ aL , - 
V\ f+ p w v%^[(ctf+pmo)^ + eov{F,P)6r V ]V?^ 



2 

Let us introduce new variational functions in the coordinate representation 



i / dr[V<p(r,P)] 2 . (65) 



In order to clarify the physical meaning of these functions let's note that if the internal degrees of freedom described 
by the bispinors Uq* sv ; Vq* av were not taken into account, the wave function P) could be represented as 

*(f,P) = e^%[(r),P = 0)]. 

However, in the problem under consideration we can't use this representation directly because of the mixing up of 
translational and spin variables. 

So, the variation of the functional (|65|l by the classical field leads to 

E{P) = J dr{* + (r,P)[(aP-iaV + /3m ) + e ^(r,P)]*(f,P)- 
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* +c (t% P)[(-aP - iaV + f3m ) + e ~ip(r, P)]* c (r, P); 

p(r,P) = e / _fL_[4+(^,P)^,P)-^+=(^,P)^,P)]; 
7 |r-r| 

y dr[*+(r,P)*(r,P) + * +c (r',P)1' c (r',P)] = 1. (67) 

As it was mentioned above, there is quite a close analogy between the nonperturbative description of QED and the 
theory of strong coupling in the "polaron" problem. Therefore it is of interest to compare the obtained functional 
116 71 with the results of various methods of taking into account the translational motion in the "polaron" problem. 
The simplest one was used by Landau and Pekar [32] who introduced the Lagrange multipliers in the form 



J(P) = J(P = 0) + (PV), (68) 

with the functional J(P = 0) referring to the resting "polaron" and the Lagrange multiplier Vi having the meaning 
of the components of the quasi-particle average velocity. 

We can see that the obtained functional l|67|) has the same form as (|68|l if the relativistic velocity of the excitation is 
determined by the formula 



V = J dr{*+(r,P)(a)*(r,P) + * +c (r,P)(a)* c (r,P)}. (69) 
It corresponds to the well known interpretation of Dirac matrixes Q. 

Varying this functional by taking into account the normalized conditions leads to the following equations for the 
wave functions 



{(aP - iaV + j3m ) + e <p(f, P) ~ £(P)}*(?% P) = 0; 

{{aP + iaV - /3m ) - e ^(r, P) - E(P)}^ c {f, P) = 0. (70) 

These equations show that distinct from the "polaron" problem |22| the translational motion of the quasi-particle 
with the momentum P in our case is related to its internal movement described by the coordinate r only by means of 
spinor variables only. The physical reason for such separation of variables is explained by the fact that in QED the 
self-localized state is formed by the scalar field, and its interaction with the particle doesn't include the momentum 
exchange. In order to find the analytical energy spectrum E(P) the system of non-linear equations J7U| should be diag- 
onalized relative to the spinor variables. The possibility of such diagonalization seems to be a non-trivial requirement 
for the nonperturbative QED under consideration. 

The solution of the equations l|70(l can be found on the basis of the states for which the dependence on the vector 
q in the wave packet amplitudes U^ s , V^ s remains the same as it was in the motionless "physical" electron. However, 
the relation between the spinor components of these functions can be changed. But as the self-consistent scalar 
potential includes the summation over all spinor components it can be supposed that the potential doesn't depend 
on the momentum for the class of states in question: 



<p{f,P) = <p(T)\ p=0 . (71) 

In the coordinate representation the transformation of the spinor components of the wave functions satisfying the 
Equations (|70|l takes place because of the dependence on the momentum. It is possible that there are some general 
arguments based on the hamiltonian symmetry that make easier the diagonalization of Equations (|70|) . However, we 
have constructed the solution by exhausting various linear combinations of the wave functions found in Section 3 for 
a resting electron. These functions correspond to the degenerated states in the case of P — but are mixed for a 
moving electron. It was found that there is the only normalized linear combination satisfying to all the necessary 
conditions of self-consistency : 

*(r, P) = L{P)^{r) + K(P)^ c (r); 
y c {r 1 P)=L 1 {P)t> c {r)+K 1 {P)^{r)- 

\L\ 2 + \K\ 2 = l^l 2 + \K X \ 2 = 1. (72) 
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The condition (|71() according to which the potential doesn't depend on the momentum for the excitation with the 
energy E(P), is fulfilled if the coefficients are related as 

Li = -K; K x = L. (73) 
These conditions are also consistent with Equations (|70|l for wave functions. 

This means that the " physical" electron moves in such a way that its states are transformed in the phase space of the 
orthogonal wave functions (|38|l . (|42() . I|46|l . (|47|l but the amplitudes of its "internal" charge distributions are not changed. 
These results, however, are valid only for the case of neglecting the interaction with the transverse electromagnetic 
field. 

Substituting the superpositions 172fl to Equations (|7U[) we use the following relations: 

(aP)W) = ( if ^S^ M ) ; [*PW{r) = ( 9 }^f^ M 
V 9{ r )(< J P) il l/2,0 : M J V -«7i(r)((7^)!!i/2 ) o,M 

H aV + p mo + e 09)m = E(0) ( ftffi^ ) 5 
H«V + fa + eo^W = -E(0) ( ) , (74) 

and similar formulas for the functions \I/(r*); \& c (r); ai are the Pauli matrixes. 

For Equations (|70(l to be fulfilled for any vector f and it is necessary to set the coefficients of spherical spinors equal 
to the same indexes I. The corresponding radial functions are proved to be the same under these conditions and the 
following system of equations for the spinors Xo l ^ s obtained 

tL(<jP)xt + K(E + E )xt = 0; iK{aP)xt + L(E - E )xt = 0; 
L(aP)xt + iK ( E + E o)Xo = 0; K(aP)xt + iL{E - E )xf = 0; 
L 1 {aP)xi - iKx{E - E ) X i = 0; K 1 {aP) X a - iL x (E + E )xo = 0; 

iLi {aP)xi — Ki(E — E )xi = 0; %K X {aP)xo ~ Li{E + E ) X o = 0. (75) 

Spin variables in Equations H75J) are also separated. In order to show it one can use, for example, the relation 
between the coefficients resulting from the 4th equation in (|75|) . 

K(aP)xt 

Al 1 



L(E-E ) 

in the first one of these equations. As the result there exists a non-trivial solution of these equations for two 
branches of the energy spectrum 



E e , p = ±yjtf* + P2 ; ( 76 ) 

referring to the electron and positron zones, respectively (Fig.l). The same expressions can be obtained for all 
conjugated pairs of the equations in J7SJ). The coefficients in the wave functions (|72|l can be found on the account of 
the normalization condition: 

"-is*- P ■ K e = -Lt Ee ~ E ° 



V^ + OEe-^o) 2 ' VP 2 + We - Eo) 2 ' 

LP = K\ = P =; K* = -L\ = ; Ee + E ° = . (77) 

VP 2 + (E e + Eo) 2 ^/P 2 + [Ee + E ) 2 

One can see that this set of coefficients coincides with the set of spinor components for solving the Dirac equation 
for a free electron with the observed mass m = Eq. Thus, the results of this section show that the "internal" structure 
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of the " physical" electron (positron) considered in this paper is consistent with the experimental energy dispersion 
(|61|l for a real free particle due to the relativistic invariance of the Dirac equation. 

However, in order to make the interpretation of the wave packet (IB2I) as the state vector for a "physical" particle 
complete we should consider one more problem that has not been discussed yet. The matter is that the exact vectors 
of the one-particle states of EPF should not only be normalized according to ljf)5|) but be also orthogonal 



<$i(A)|«i(P)>=«(P-A). (78) 

It is evident that the approximate representation (|62(l for the one-particle state vector doesn't satisfy this condition 
. The reason for non-orthogonality is that the variational principle J7UJ)used for determining the energy spectrum 
of quasi-particle excitations allows to satisfy the exact law of the momentum conservation only on the average as 
in Equation JHSJ- At the same time the orthogonality condition (|78|l holds true for the exact eigenfunction of the 
momentum operator. 

The rigorous method of taking into account the translational symmetry in the theory of strong coupling for the 
"polaron" problem was elaborated in the works of Bogolubov and Gross (2^. Let's remember that this method 

was based on the introduction of the collective variable R conjugated to the total momentum operator P, the canonical 
character of the transformation caused by three new variables Ri being provided by the same number of additional 
conditions imposed on the other variables of the system. In the "polaron" problem the quantum field interacting with 
the particle contributes to the total momentum of the system. It allows to impose these conditions on the canonical 
field variables [25| and the concrete form of the variable transformation is based essentially on the permutation 
relation for the operator of the boson field . 

The considered problem has some specific features in comparison with the "polaron" problem. Firstly, the formation 
of the one-particle wave packet is the many-particle effect because this packet includes all initial states of EPF as the 
fermionic field. Secondly, its self-localization is provided by the polarization potential of the scalar field which doesn't 
contribute to the total momentum of the system. Therefore we use another approach in order to select the collective 
coordinate R. Let us return to the configuration representation in the hamiltonian l|19l) , where QED is considered 
to be the totality of N (N — > oo) point electrons interacting with the quantum EMF in the Coulomb gauge 



H = ^{a a [p a + e Q A(f a )} + f3 a m + e <p(f a )} 

a=l 

X - j dr(V<p(r)) 2 +Y j uj{k)h ix , 

kX 



u{k) = k\ n^=c+ x c %x ; A =1,2; 
2(f) = y - 7 ^e< A > [cejV 6, + cf e~ i%f ] (79) 

kX 



Here fi is the normalized volume; ct (cg A ) are the operators of the creation (annihilation) of quanta of the transversal 



feA 

electromagnetic field, the quantum having the wave vector k, polarization e^) and energy uj{k) = k. The sign of the 
interaction operators differs from the standard one because the parameter eo is introduced as having a positive value. 

As it was stated above the zero approximation of nonperturbative QED is defined only by a strong interaction of 
electrons with the scalar field, the interaction with the transversal field is taken into account in the framework the 
standard perturbation theory when the conservation of the total momentum is provided automatically 0. Therefore, 
while describing the quasi-particle excitation we should consider the conservation of the total momentum only for 
the system of electrons. In the considered representation it is defined by the sum of the momentum operators of 
individual particles 

N 
a=l 

P\<f>l(P) >=P|$l(P) > . (80) 

It means that in the configuration space the variable R conjugated to the total momentum is simply a coordinate 
of the center of mass and the desired transformation to new variables is as follows: 
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^ N N 

r a = R + p a : R=j-f^2r a ; p a = 0; 



a—1 a—1 
N N 



#a = -^V a = -^V R +# ; P = -iV fl ; pt--*V, a + -^^V P6 ; = 0. (81) 

b=l a=l 



The hamiltonian i|79|) with new variables has the following form 

JV 

H = ^{M-^R +A + eoApa + R)] + Pam a + e o 0(p a + R)} 



l - J dr(V0(r)) 2 + X>(*)fij£v (82) 

k\ 



It is essential to note that the matrix elements of an arbitrary operator in a new configuration representation should 
be calculated in accordance with the following norm (we introduce a special notation for that norm) 

« $x|M|$ 2 »= f dR]]dp a $l(R,{p a })M<I> 2 (R,{p a }). (83) 

J a 

Let us denote by Hq that part of the operator l|82|) which doesn't depend on the transversal EMF and describes 
the internal structure of "physical" particles . Actually the operator Hq doesn't depend on the coordinate R either, 
because of its commutativity with the operator of the total momentum of the system of electrons. This also follows 
from the well known result [lj] that in the Coulomb gauge the scalar potential could be excluded from the hamiltonian. 
As the result the operator Hq depends only on the vector differences (r a — f%) = (p a — p ) not changing with the 
simultaneous translation of all the coordinates. As the consequence, the eigenfunctions of the hamiltonian Hq depend 
on the coordinate R in the same way as for a free particle: 

N i i r 

H -> H (P) = ^{^[-P + gJ + #> TO ° + e o^(Pa)} - 2 / ^W) J ( 84 ) 

a=l J 

Further calculations consist in returning to the field representation by the variables p a in the limit (N — > oo) and in 
using the approximate trial wave packet l2"2"|) for the state vector |$i(P) > . As the result one obtains the functional 
for the energy E(P) which coincides completely with the expression (|67J) . It should be noted that the averaging of the 
operator H$(P) in the secondary quantization was calculated using the relation which follows from the well known 
density of states when EPF is quantized in the normalized volume Q 

lim — = 1. 

P 

Thus, the orthogonalized and normalized set of states for the EPF one-particle excitations corresponding to the 
"physical" electrons (positrons) is defined as follows 

|4>i O) 0P) >^ j2^m GlM J dq{U ?s (P)a% + V q ^P)btJ\0-,0;^ >; (85) 

with the norm (|83|) and the coefficient functions U^ S (P); V^ S (P) which are combined by the Fourier transformation 
of the wave functions *f?(r, P); \E' c (r, P) defined by the formulas 
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V. PERTURBATION THEORY FOR QED WITH A LARGE VALUE OF THE INITIAL COUPLING 

CONSTANT 



It is well known that the renormalizability is one of the most important features of QED and it is confirmed by 
the coincidence of its results with the experimental data. Actually it means that the calculated characteristics of the 
real electromagnetic processes don't depend on the initial parameters eo,ni but only on the observed values of e 
and m. So, we should clearly show that the effects of the nonperturbative QED refer to the internal structure of the 
"physical" electron (positron) but its interaction with a transversal electromagnetic field is defined by the standard 
perturbation theory with the fine structure constant a <C 1. As for a scalar field, this requirement is met by the 
normalized condition l|24l) . 

One more problem is to demonstrate that the integrals which define the corrections to the zero approximation of the 
nonperturbative QED are converged without introducing any additional regularizating parameters. This condition 
was one of the reason for the present work. 

To solve these problems the form of the perturbation theory which uses the basis of states corresponding to 
" physical" electron (positron) should be considered. Taking into account the results of the last section let us represent 
the QED hamiltonian (|19|) in the following form : 

H = H + Hi; 

H = J dr{: 4>*(r)[-ia(f $ + v>) + /3m ]^(r) : +e (^(r) : p(r) : -i(V£(r)) 2 } + ^ u(k)h Ux ; 

k\ 

= e J df : tfj* -(r)[aA(r '+ R)]i>{r) : . (86) 

The field operators i/j(r) and the density operator p{r) are defined by expressions lfT$ )l .l|2() [l and we the symbol : A : 
of the normal ordering is used to exclude the vacuum energy l|28f) in the EPF excitation spectrum. The fact that the 
hamiltonian depends on the coordinate R canonically conjugated to the total EPF momentum determines the local 
violation of symmetry of the system. This situation was discussed by Bogolubov |2l| for a lot of concrete physical 
systems. However, the global symmetry of the system |21| is reconstructed when the observed values are averaged by 
the coordinate R in accordance with the norm 

As distinct from the canonical perturbation theory we can't find exactly the total spectrum of eigenvectors of the 
operator Ho for the considered case of the electron strong coupling with a scalar field . In the previous sections only 
the spectrum of the one-particle excitations was found approximately in the main order of the expansion in terms 
of ag 1 power. In particular, even the calculation of the two-particle excitation spectrum requires the analys of the 
system analogous to the " bipolaron" problem in the Pekar model |23] ■ Such analysis is out of the framework of the 
present paper. Therefore we consider only such matrix elements of the perturbation operator Hi which are due to 
the transitions between the states corresponding to the free "physical" electron (positron) and an arbitrary number 
of quanta of the transversal electromagnetic field. Let us write the obvious form of these states using definition 185|) 
for the one-particle wave packet 



N u A >^ jTT^e 1 ™ / dq{U q ~ s (P)a+ + V 9lU {P)b+}-^=^\0; 0; <p(f) > . (87) 
[Z7r > J \ N kJ 



In order to be definite, let us calculate the matrix element of the operator Hj for the transition of the "physical" 
electron between the states with the momenta P and P\ together with emitting one quantum of the transversal EMF. 
The norm (|8*3j) is used when calculating this element 



T A (P, k) =« $i(Pi), fc, A|#/|$i(P), »= S(P - P 1 - K)~/ X (P, k); 



7 (P,fc) 



'2kQ. 



EE 



df 



dq 



dq 1 



;{UlAPl)Ug,s(P)u^ s ,^de X )^U q ^ 



(2tt) 3 / 2 J (2tt) 3 / 2 l q ' s K ' q, y ' g s ^ 

^ )S ,(A)% s (P)^ sV («e A ) / ,,^- fs Je i(? - 9 "'-" )f? - 



(88) 



In order to calculate the renormalized vertex function 7(P, k) let us find the under integral expression through the 
spinor wave functions of one-particle states in the coordinate representation by means of the relations l|66|l 
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j X (P,k) = / df{**(f,A)(ae A )*(f,P) - * c '*(r, A)(ae A )* c (f, P)}e- i%p \ 



P x =P-k. (89) 



Let's substitute the linear combinations (|72|l in this expression and use the analytical form of the functions 
^, vp c for the resting "physical" electron given in Eqs. (|2HJl , (El , (EH i • Then the transition to the two- 
component spinors should be written as follows: 

7A (P,fc) = i-^g= | dr{[ 5 2 (r)( Xo +*(^ A )x + ) - / 2 W(xr (^) X f)][i^(Pi)i(P) - 

L*(P 1 )X(F)]-[.9?(r)(xr*('?e A )xr)- 
f?(r)(Xo*(^ X )Xo)][K!(Pi)Li(P) - Lt(P 1 )X 1 (P)]}e-* (90) 

In accordance with the results of Section 3 the radial wave functions of the internal structure of the "physical" 
electron are localized in the domain r < (6mo) _1 - It means that if the vertex function is calculated with the photon 
momentum less than a certain finite but large enough value 

k < ko ~ bm Q ~ 328m ~ 167(MeV), (91) 

the exponent in the integrals H9Q(I can be substituted by unit exp(— ikr) ~ 1. Besides, the spinor and angular 
variables can be excluded by means of Eq. (|72(l for the spinor wave functions describing the internal state of the 
"physical" electron. Thus, for example, 

K(P)(^) X + = -^ { aP)(ae>) X + = ^- ( /% A ) Xo +. 

Taking into account the normalization of the spinors xji~ one can find that the vertex function has the following 
form when the condition (|91|) is fulfilled 



V2W, E(P)+E P(Pi) + P 

dr{[9 2 (r) + f 2 (r)]-[g 2 1 (r) + fUr)}} = 

= L(P)L(A)[-S^ , ( ? gA) ]• (92) 



V2M7 J B(P) + m P(Pi) 

We have also used the normalization conditions for the radial functions and the relations l|51|l which define the 
observed charge e and mass m of the "physical" electron. The obtained expression for the vertex function completely 
coincides with the result of the standard perturbation theory [Llj . This result is not obvious a priori and it seems to 
be a serious argument in favour of the nonperturbative QED. It also shows that after the renormaliation our approach 
leads mainly to the same results for the observed characteristics of the electromagnetic processes as the standard form 
of the perturbation theory in QED. 

And finally, the last question we are are going to consider briefly in the present paper is concerned with the 
regularization of the diverged integrals in QED. It is evident from expression (|90|) that the behavior of the renormalized 
vertex function 7(P, k) at large k is defined by the Fourier-image of the radial functions connected with the form-factor 
of the charge distribution in the "physical" electron. If the variational solutions (|53|l are used, then the asymptotic 
behavior of the vertex function at large k 3> ko is defined by the following dependence 
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and converges rather quickly to zero . 

Thus, the cut off momentum k max ~ k Q appears naturally in the nonperturbative QED and it has quite a def- 
inite physical meaning. It is well known (for example, Q) that the introduction of such momentum provides the 
regulariztion of all integrals in QED. However, it should be stressed that the wave functions describing the charge 
distribution in the "physical" electron are not completely spherically symmetrical . Therefore, the effective cut off 
of the integration domain in the integrals over the momenta could be not spherically symmetrical either unlike of 
the usual methods of regularization. But these corrections should be rather small because of a quick decrease of the 
functions under the integrals at the domain boundary in accordance with l|93|l . The concrete estimation of these 
corrections is connected with rather cumbersome calculations and is out of the framework of the present work. 



VI. CONCLUSIONS 



Thus, a new logically consistent solution of quantum electrodynamics equations is constructed in the present paper. 
It is based on the supposition that the existing form of QED is a closed and self-consistent theory of the interacting 
electron-positron and electromagnetic fields. 

The main advantage of the obtained solution is the possibility to connect in the finite form the observed and initial 
parameters of QED without introducing any additional parameters. This solution, of course, doesn't exclude the 
interaction of EPF and EMF with other quantum fields, but demonstrates the possibility of the finite renormalization 
in QED itself without taking into account these fields and other external factors. 

According to our results the non-renormalized form of QED corresponds to the theory of strong interaction between 
electron-positron and scalar electromagnetic fields. The one-particle excitation of EPF leads to a local modification 
of the vacuum states of the system, with quite a strong classical self-consistent field appearing as the result. In the 
field of this potential the "physical" electron (positron) is formed as the quasi-particle excitation consisting of two 
strongly coupled charge distributions of opposite signs. The effective masses of both distributions are large. But the 
integral charges of both distributions are almost equal to each other and their binding energy is comparable with the 
sum of their masses. As the result of the compensation of these large values the observed values for the charge and 
mass of the "physical" electron (positron) appear. It is shown that the electromagnetic processes for the "physical" 
electron-positron field are defined by the renormalized charge and mass of the electron. But the internal structure of 
the "physical" electron defines the unobserved cut off momentum in the vertex function. It is also shown that the 
states and energy spectrum of the "physical" particles satisfy the conditions of the relativistic invariance. 

Certainly, the calculations given in our paper are approximate and should be developed in various directions. In 
particular, it is important to formulate the nonpertubative QED in a covariant form on the basis of the Lorentz 
gauge, to find a decisive evidence that the results of nonperturbative and standard QED are equivalent after the 
renormalization, to investigate the effect of other fields on the obtained solution, to consider the spectrum of two- 
particle excitations. The formulation of the nonperturbative QED on the basis of the path integrals is of special 
interest. It is impossible to analyze all these problems in the framework of one paper but it seems to us that the physical 
picture of the renormalzation considered here could be the foundation for further discussions and investigations in 
the above mentioned directions. 
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CAPTION FOR FIGURES 
FIG.l 

Scheme of the one-particle levels for the vacuum (£(°)(P)) and "physical" (E(P)) electron-positron field. 



